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Abstract 

We report measurements of the D^, D^, and Df meson lifetimes using 

3.7 fb^^ of e~^e~ annihilation data collected near the T{AS) resonance with 

the CLEO detector. The lifetimes of the D^, D^ , and D^ mesons are mea- 

^ ■ sured to be 408.5 ± 4.lt|4 fs, 1033.6 ± 22.llJ-j fs, and 486.3 ± IS.Otg;? fs, 

respectively. The precisions of the charm meson lifetimes reported here are 

— ^ . comparable to those of the best previous measurements, and the systematic 

fSJ \ errors are very different, br a single experiment we find that the ratio of the 

^P ' Dt and D^ lifetimes differs from one by more than 4.5 standard deviations. 
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The systematics of charm hadron hfetimes have played a central role in understanding 
heavy quark decays |[l|]. In this Letter we report new measurements of the lifetimes of the 
D^, D~^, and D^ mesons. These charm meson ground states differ in the identity of the light 
antiquark, i.e., the D^, D", and D^ mesons are cd, cu, and cs states, respectively. Although 
the weak decay of the charm quark is responsible for the decays of all three charm mesons, 
differences in the lifetimes indicate that the identity of the light antiquark also influences the 
rates of decay. The large ratio of the D~^ and D^ lifetimes {t£)+/tdo ~ 2.5) arises primarily 
from destructive interference between different quark diagrams that contribute significantly 
only to decay of the D^ |jl|. This interference, as well as a number of smaller effects, which 
can cause the D'^ and D^ lifetimes to differ, appear in a systematic expansion, in inverse 
powers of the charm quark mass, of the QCD contributions to the charm decay amplitudes [|l|. 
The results described in this Letter indicate that the ratio of the Df and D^ lifetimes differs 
significantly from one, providing a quantitative challenge for the theoretical description of 
charm meson decays. These data were obtained in an e+e^ colliding beam environment, 
where the event topologies and backgrounds are very different from those encountered in 
the high energy fixed target experiments ||^ that have recently provided the most precise 
measurements of charm hadron lifetimes 0. 

The results described in this Letter are based on an integrated luminosity of 3.7 fb~^ of 
e~^e~ annihilation data recorded with the CLEO II. V detector near the T[4S) resonance at 
the Cornell Electron Storage Ring (CESR). This luminosity corresponds to approximately 4.4 
million recorded e'^e" —>■ cc events. The CLEO II detector has been described elsewhere [Q. 
The major component of the CLEO II. V upgrade completed in November 1995 is the SVX, 
the first multi-layer silicon vertex detector operating near the T(4S') energy ||^. The SVX 
consists of three concentric layers of 300 /im thick, double-sided silicon strip detectors to 
measure the xy and rz coordinates ^ of charged particles. The three layers are at radii of 
2.35, 3.25, and 4.75 cm. There is a total of 0.016 radiation lengths in the material in the SVX 
and the beryllium beam pipe whose inner radius is 1.875 cm. For this detector, the average 
"signal-to- noise" ratio for charged particles at minimum ionization is 15:1 for the xy view 
and 10:1 for the rz view and the efficiency to have two or more SVX hits simultaneously in 
both views is 95% per track. The impact parameter resolutions as functions of momentum 
p (GeV/c) are measured from data to be a^y = 19 © 39/(psin^'^ 6*) yum and (at 6 = 90°) 
arz = 50 © 45 /p fim |^. The Monte Carlo simulation (MC) of the CLEO detector response 
is based upon GEANT 0. Simulated events are processed in a similar manner as the data. 

We reconstruct D mesons in the decay modes D° — *> K'n'^, D^ -^ K'tt^tt^, D^ — >■ 
K-'K+'K-'K+, D+ -^ K-7r+7r+, and D+ -^ (j)7r+ with (p -> K+K'. In this Letter, "D" refers 
to D°, D^, and Df mesons and reference to the charge conjugate state is implicit. The 
charged D daughters are required to have well reconstructed tracks and to have particle 
identification information from specific ionization {dE/dx) and time-of-fiight consistent with 
the D daughter hypothesis. Charged tracks forming a D candidate are required to originate 
from a common vertex. Neutral pions are reconstructed from photon pairs detected in the 
electromagnetic calorimeter. The photons are required to have an energy of at least 30 
(50) MeV in the barrel (endcap) region and their invariant mass is required to be within 



three standard deviations of the nominal vr" mass. The vr*^ momentum for D^ -^ K tt'^tt 







is required to be greater than 100 MeV/c. For background suppression, a soft pion vr^ 
(vr^) is required to form a D*~^ with the D candidate for the D^ {D~^) decay modes. The 



£)*+ _ £)0 (^/}+) mass difference is required to be within 800 (1400) keV/c^ of the nominal 
value 0. No such requirement is made for the decay D^ -^ 07r+, where the requirement 
that the K^K~ invariant mass be within 6 MeV/c^ of the mass substantially reduces the 
background contribution. In the decay D^ -^ 'P'^^ ■, the angle between one of the kaons and 
the pion in the rest frame of the meson follows a cos^ 9 distribution. Since the combinatorial 
background for this decay is distributed uniformly, we require | cos^| > 0.4. The D*^ and 
the D^ momenta are required to be greater than 2.5 GeV/c. The mass distributions for the 
D candidates (after subtracting the nominal D mass values p[) are shown in Fig. ||. The 
numbers of reconstructed D mesons N£,, given in the figures, are obtained from fits of the 
mass distributions to two Gaussians over a linear background. The background fractions in 
the mass regions within ± 16 MeV of the nominal D mass values are 1.2% [D^ — > ir~7r+), 
4.9% (DO ^ ir-vr+TrO), 10.0% (D^ ^ R- ti+ ti' ti+) , 12.2% (D+), and 13.8% (D+). 
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M(D)-Mp(GeV/c^) 

FIG. 1. Masses of charmed meson candidates M{D) minus the nominal masses M^, for (a) 
D^ -^ i^~7r+, (b) D° -^ K-TT+7r°, (c) D° -^ i^-7r+7r-7r+, (d) D+ -^ A'^7r+7r+, and (e) D+ -^ 07r+. 
The data (solid squares) are overlaid with the fit to two Gaussians with the same mean over a linear 
background (solid line). The fitted background is indicated by the dashed line. 



The dimensions of the CESR luminous region (beam spot) are known from the machine 
optics to be about 1 cm along the beam direction (z), 7 /xm in the y direction, and about 
350 /xm in the x direction. The centroid of the beam spot is determined for each CESR 
fill. The charm hadrons are produced approximately back-to-back at the primary vertex 
(interaction point). In the laboratory frame the selected D^, D~^, and D^ mesons have an 
average momentum of 3.2 GeV/c and average decay lengths of 200 fim, 500 /xm, and 240 /xm. 



respectively. 

The decay vertex r^ and the momentum vector p^^ of each D meson candidate are 
reconstructed in the xy plane. The decay vertex resolution along the D flight direction is 
80 — 100 /zm depending on the decay mode. The interaction point rjp is reconstructed by 
extrapolating the D momentum back from the decay vertex to the beam spot. We calculate 
the projected decay length /dec from the displacement in the xy plane from rjp to the D 
decay vertex, /dec = (i"d — rip) ■ p^. We then calculate the proper time of the D meson decay 
from t = mE>ldec/cPD using the appropriate PDG [0 world average for mo- The proper-time 
distributions for the D candidates are shown in Fig. ^j. 

The D meson lifetimes are extracted from the proper-time distributions with an unbinned 
likelihood fit. The likelihood function is 

^\'Di Jhgi 'bg) *^) /mis) ''^mis) /wide J 






Psig,rEit'\Tn) + (1 - p.ig,.) [/bg^(tlrbg) + (1 - /bg)(5(t')] 

"^ V ' "^ V ■■ 

signal fraction background fraction 



,1 ~ /mis ~ fwidc)G{ti — t |oO"t,j) + fmisGyti — t IcTjnis) + fwidcG{ti — t IcTwide) 



proper— time resolution mismeasured fraction 

where the product is over the D meson candidates, G(t\a) = exp(— t^/2(T^)/\/27ro", and 
E{t\T) = exp{—t/T)/T. We fit the proper-time distributions for the different decay modes 
separately. In these fits, each D meson candidate is assigned a signal probabihty psig,i 
based on its mass. The signal probabilities are derived from the (independent) fits of the D 
mass distributions to the sum of two Gaussians with the same mean and a linear background 
function. The seven parameters of the lifetime fit are Tn, /bg, Tbg, S, /mis, amis and /„ide- The 
parameter td is the D meson lifetime. The background proper-time distribution is modeled 
by a fraction /bg with a background lifetime Tbg and a fraction with zero lifetime. In order to 
estimate the background properties, we fit the candidates in a wide region of ±40 MeV/c^ 
around the nominal D mass. Each candidate is weighted in the fit according to its proper- 
time uncertainty at^i- The fit allows for a global scale factor S that modifies the calculated 
proper-time uncertainty. The fits yield S* ~ 1.1 for all modes. For a small fraction of 
mismeasured candidates /mis, the fitted uncertainty Sat^i underestimates the true uncertainty. 
This is a result of track reconstruction errors such as hard multiple scattering or the use of 
an SVX noise hit in the track fit. In the fit, we account for the mismeasured candidates 
with two Gaussians. The fit parameters associated with the mismeasured candidates are 
the fraction of events in each of the Gaussians /mis and /„ide and the width of one of the 
Gaussians o"mis- The width of the other Gaussian (uwide = 8 ps) is fixed. The results of the 
unbinned likelihood fits are superimposed on the proper-time distributions shown in Fig. ^ 
From the fits we obtain tdo = 411.1±5.7 fs iK-n+), 395.2±8.1 fs (K-7r+7r°), 416.3±8.6 fs 
(K^Tf^TT^Ti^), T£,+ = 1033.6 ± 22.1 fs, and r^+ = 486.3 ± 15.0 fs, where the uncertainties 
are statistical only. The correlation coefficients between the D lifetime and the other fit 
parameters are typically near 0.1, and the largest is 0.28. All of these fit results have been 
corrected for small biases observed in the measurements of the D lifetimes in simulated events 

6 
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t(ps) 
FIG. 2. Proper-time distributions of charm meson candidates within ±16 MeV/c^ of the nom- 
inal D mass for (a) L>° -^ K-7r+, (b) D° -^ K-ir+ir^, (c) D^ -^ K-tt+t^-tt^, (d) Z)+ -^ K-tt+tt+, 
and (e) Df -^ </)7r+. The data (soHd squares) are overlaid with the result from the unbinned likeli- 
hood lifetime fit (solid line) . The proper-time spectra of the background candidates obtained from 
the fits are indicated by the shaded area. 



of -3.0 ± 0.9 fs {K-n+), 2.4 ± 2.3 fs [K-tt+tc^), -2.0 ± 2.2 fs [R- 71+71-71+), -2.9 ± 6.6 fs 
{D+), and -0.6 ± 2.4 fs {D+). The D^ lifetime tdo = 408.5 ± 4.1 fs (combined) is obtained 
as the weighted average of the three measurements using statistical uncertainties only. The 
individual measurements are weighted by their inverse relative uncertainty (r/a,-)^ |]10| . 

The large samples of reconstructed charm mesons permit a number of consistency checks, 
including varying the D candidate mass region, measurement of the background properties 
in the D mass sidebands, and division of the data samples in several key variables such as 
azimuthal angle, polar angle, momentum of the D candidate, and data taking period. No 
statistically significant effect is found in any of these variables. 

The systematic uncertainties for the D meson lifetimes are listed in Table |I| and are 
described below. They can be grouped into three categories: 

Reconstruction of the D decay length and proper time. Errors in the measurement of the 
reconstructed decay length can be due to errors in the measurement of the decay vertex, 
the global detector scale, and the beam spot. The bias in the decay vertex position is 
estimated to be (0.0 ± 0.9 /im) from a "zero-lifetime" sample of 77 — »• 7r"'"7r~7r"'"7r~ events. 
This corresponds to a measured proper-time uncertainty of ±1.8 fs. In addition, the vertex 
reconstruction is checked with events with interactions in the beam pipe with a relative 



TABLE I. Systematic uncertainties for the D meson lifetimes in fs. The systematic uncertainties 
for the three D^ modes are weighted with the same weights as the fitted D^ hfetimes. 



Uncertainty 




Z)0 




D^ 


D+ 


Dt 




K-TT+ 


K--K+-K^ 


K vr+vr 7r+ 


combined 


K"7r+7r+ 


(/)7r+ 


Decay vertex 


±2.0 


±2.0 


±2.0 


±2.0 


±2.8 


±2.1 


Global detector scale 


±0.1 


±0.1 


±0.1 


±0.1 


±0.1 


±0.1 


Beam spot 


+0.3 
-0.1 


+2.1 
-0.0 


+0.3 
-0.2 


+0.8 
-0.1 


+1.3 
-1.1 


+0.7 
-1.1 


D meson mass 


±0.1 


±0.1 


±0.1 


±0.1 


±0.3 


±0.1 


D meson momentum 


+0.2 
-0.0 


+0.1 
-0.2 


+0.3 
-0.1 


+0.2 
-0.1 


+0.6 
-0.0 


±0.1 


Signal probability 


+0.4 
-0.1 


+0.1 
-0.2 


+0.1 
-0.2 


+0.3 
-0.1 


+1.2 
-8.1 


+1.3 
-1.8 


t - M{D) correlation 


±0.6 


±0.6 


±1.0 


±0.7 


±1.7 


±1.5 


Large proper times 


±1.2 


±3.4 


±0.2 


±1.5 


±0.3 


±0.5 


Background 


±0.5 


±2.4 


±3.0 


±1.5 


±6.3 


±2.9 


MC statistics 


±0.9 


±2.3 


±2.2 


±1.6 


±6.6 


±2.4 


Total 


+2.7 
-2.6 


+5.t) 
-5.2 


±4.4 


+3.5 
-3.4 


+ 9.9 
-12.7 


+4.9 
-5.1 



uncertainty of ±0.2%. The sums of these uncertainties in quadrature yield the systematic 
uncertainties due to the decay vertex measurement. The global detector scale is measured to 
a precision of ±0.1% in surveys and confirmed in the study of events with interactions in the 
beam pipe. The changes in the lifetimes due to the variation (±2 /im) in the vertical beam 
spot position and height are another source of systematic error, since the interaction point 
is calculated from the beam spot and the reconstructed D momentum and decay vertex. 
Statistical uncertainties for the D masses [Q and the D momentum measurements lead to 
systematic errors since these quantities are used to convert the decay length into proper 
time. 

Lifetime fit procedure. This category includes uncertainties in the candidate signal prob- 
abilities, the impact of candidates with large proper times, the correlation between proper 
time and D meson mass, and the proper-time properties of the background. The signal prob- 
ability assigned to each candidate in the lifetime fit has a statistical uncertainty, and these 
statistical uncertainties lead to systematic uncertainties in the fitted lifetimes. We estimate 
these systematic uncertainties by coherently varying the signal probability of each candidate 
by its statistical uncertainty and repeating the fits. A correlation between the measurements 
of the proper time t and the D candidate mass M{D) can be a source of systematic uncer- 
tainty. We measure this correlation in simulated events to estimate the associated systematic 
uncertainty. Charm meson candidates with large proper times are an additional source of 
systematic uncertainty. These candidates are modeled by the wide Gaussian in the proper- 
time fit. Alternatively, the wide Gaussian component is omitted from the likelihood function 
and candidates in a restricted proper-time interval are fitted. The systematic uncertainties 
due to candidates with large proper times are estimated from the variations of td with the 
width of the wide Gaussian and the differences in the results between the fits with different 
proper-time intervals. This systematic uncertainty is small for decay modes with three or 



more charged D daughters for which the requirement of a well-reconstructed vertex greatly 
reduces mismeasurements. We estimate the systematic uncertainty due to charm and other 
backgrounds that might populate the D mass peaks differently than they populate the D 
mass sidebands, 20 MeV/c^ < \M{D) — Md\ < 60 MeV/c^. Some possible sources of such 
backgrounds are a background in the D^ sample from D^ ^ K^7i^tt~ decays where one vr" 
is misidentified as a K~ , and backgrounds from D^^^^ decays in the D^^'^^ sample caused by 
adding or missing a charged pion. 

Checking the algorithms with simulated events. Charm meson candidate selection re- 
quirements can cause systematic biases in the lifetime measurements. We estimate these 
biases with simulated events and correct for the biases as described above. We include the 
statistical uncertainties in the measured lifetimes from the samples of simulated events as 
systematic uncertainties in the results. 

The total systematic uncertainties in the D^ lifetime measurement are obtained by com- 
bining the contributions from the three reconstructed D^ decay modes. The contributions 
from the decay length measurement and the detector size are assumed to be completely 
correlated and all other contributions are assumed to be uncorrelated. The total systematic 
uncertainties are obtained by adding the individual contributions in quadrature. 

In summary, we report measurements of charm meson lifetimes from 3.7 fb~^ of in- 
tegrated luminosity recorded with the CLEO detector. The measured D lifetimes are 
tdo = 408.5 ± 4.1+11 fs, td+ = 1033.6 ± 22.lli|^ fs, and Tjy+ = 486.3 ± 15.0l|? fs, 
where the first uncertainties are statistical and the second systematic. These results im- 
ply t^+/t£)0 = 1.19 ± 0.04, a difference of more than 4.5 standard deviations in a single 
experiment. The charm meson lifetimes reported in this Letter are comparable in precision 
with the best previous measurements 0, and the systematic errors are very different. 
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